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ABSTRACT 


Synthetic  aperture  sonar  beamforming  and  signal  proc¬ 
essing  relies  on  properly  steering  and  focusing  the  aper¬ 
ture  beam  pattern  in  order  to  co-phase  all  the  received 
signals.  Due  to  the  effects  of  motion  in  the  synthetic 
aperture  sonar  problem,  the  propagation  path  between  the 
transmitter,  discrete  point  scatterer,  and  the  receiver  is 
time  varying.  Traditionally,  simple  approximations  are 
used  to  determine  these  propagation  ranges  and  angles  of 
incidence  and  scatter.  Methods  to  determine  these  ranges 
and  angles  exactly  may  significantly  improve  array  gain 
and,  therefore,  target  detection. 

This  thesis  investigates  improvements  to  SAS  signal 
processing  algorithms  using  exact  methods  for  the  calcula¬ 
tion  of  the  time-varying  ranges  between  transmitter  and 
discrete  point  scatter,  and  between  discrete  point  scatter 
and  receiver,  and  the  phase  angle  of  the  scattered  acoustic 
signal  incident  upon  the  receiver.  Using  computer  simula¬ 
tions,  exact  range  and  angle  calculations  were  performed 
for  different  scenarios  and  compared  to  ranges  and  angles 
determined  using  standard  approximations.  The  exact  ranges 
were  then  used  to  determine  incident  phase,  and  were  again 
compared  to  the  approximate  methods.  Comparison  of  the 
exact  and  approximate  methods  was  based  on  range  estimation 
error  and  percentage  error.  Improvements  in  synthetic 
aperture  array  gain  using  exact  phase  weights  based  on 
exact,  time-varying  range  solutions  are  proposed. 
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EXECUTIVE  SUMMARY 


The  goal  of  this  thesis  was  to  demonstrate  improve¬ 
ments  to  Synthetic  Aperture  Sonar  (SAS)  beamforming  algo¬ 
rithms.  SAS  systems  are  currently  in  development  to  sup¬ 
port  the  Navy' s  Long-Term  Mine  Reconnaissance  System 
( LMRS ) .  The  LMRS  will  use  a  submarine-launched  autonomous 
underwater  vehicle  to  create  a  synthetic  array  for  high- 
resolution  mine  imaging  and  classification. 

Currently,  SAS  signal  processing  techniques  make  vari¬ 
ous  assumptions  and  simplifications  in  the  determination  of 
the  signal  propagation  ranges  and  phase  shift  corrections 
required  for  accurate  acoustic  modeling.  These  approximate 
techniques  simplify  the  steering  of  the  beam  pattern  of  the 
sonar  aperture  but  also  unnecessarily  introduce  errors  that 
affect  image  quality. 

Beamforming  techniques  have  been  developed  that  can 
exactly  determine  the  parameters  required  to  steer  the  beam 
pattern  of  a  sonar  array.  Using  this  exact  method,  acous¬ 
tic  signal  propagation  in  the  ocean  for  a  SAS  system  was 
simulated.  This  simulation  was  conducted  using  stationary 
targets  to  replicate  the  problem  of  water-borne  mine  detec¬ 
tion  and  classification.  Water-borne  mines  include  buoy¬ 
ant,  tethered,  and  bottom  ocean  mines.  The  primary  focus 
of  these  simulations  was  to  demonstrate  the  improvement  in 
SAS  beamforming  that  can  be  realized  through  exact  propaga¬ 
tion  range  calculation.  With  exact  range  estimation,  exact 
phase  weights  are  determined  and  used  to  improve  array 


xiii 


gain . 


An  analysis  of  the  limitations  of  the  most  common  ap¬ 
proximations  is  presented.  The  exact  beamforming  methods 
have  no  restrictions  on  their  application,  and  this  is 
demonstrated  through  simulation.  The  simulation  of  the 
propagation  of  signals  from  transmitter  to  target  and  then 
to  the  receiver  allows  direct  comparison  of  the  standard 
approximate  methods  with  the  new  exact  methods.  Estimation 
of  signal  propagation  range,  travel  time,  angles  of  inci¬ 
dence  and  scatter  at  the  target  and  angles  of  incidence  at 
the  receiver  are  all  improved.  Finally,  fast  Fourier 
transform  beamforming  is  conducted  and  the  performance  of 
the  exact  and  approximate  methods  is  compared  using  a  sim¬ 
ple  signal  composed  of  three  different  sinusoids. 

This  research  demonstrates  that  simple  approximation 
techniques  traditionally  used  to  estimate  SAS  propagation 
ranges  are  not  useful  in  current  SAS  projects  that  employ 
high  frequency  signals  and  long  target  ranges.  This  re¬ 
search  also  shows  that  SAS  propagation  ranges  can  be  accu¬ 
rately  approximated  using  the  binomial  approximation  propa¬ 
gation  model  presented  in  this  thesis.  Although  the  bino¬ 
mial  approximation  propagation  model  is  accurate  under  most 
conditions,  this  model  places  some  restrictions  on  target 
range  and  platform  speed.  The  effects  of  operating  outside 
the  binomial  approximation  model  limits  are  demonstrated. 

It  is  shown  that  implementation  of  the  exact  propagation 
range  method  does  not  significantly  add  to  the  complexity 
of  the  beamforming  algorithm.  The  exact  method  introduces 
no  error  and  places  no  restrictions  on  target  location  and 
platform  speed.  Use  of  the  exact  propagation  model  is 
recommended  to  allow  maximum  accuracy  and  flexibility  as 
SAS  systems  evolve.  Further  investigation  will  determine 


the  improvements  that  can  be  realized  by  applying  the  exact 
time-varying  angles  of  incidence  and  scatter  to  the  evalua¬ 
tion  of  target  scattering  functions. 
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INTRODUCTION 


I . 


This  chapter  provides  an  overview  of  synthetic  aper¬ 
ture  sonar  concepts.  The  modeling  framework  used  to  simu¬ 
late  signal  propagation  through  the  ocean  medium  is  also 
introduced. 

A.  SYNTHETIC  APERTURE  SONAR 

Synthetic  Aperture  Sonar  (SAS)  systems  attempt  to  ex¬ 
ploit  the  benefits  that  can  be  obtained  from  very  long 
sonar  arrays  while  maintaining  the  actual  physical  size  of 
the  array  to  be  relatively  small.  This  is  an  attractive 
technique  because  it  would  allow  deployment  of  Unmanned 
Underwater  Vehicles  (UUV)  to  detect,  localize  and  classify 
mines  in  the  littoral  environment.  Vehicle  size  is  limited 
by  the  launching  capability  of  the  controlling  ship  or 
aircraft  and  a  need  to  conduct  covert  surveillance  and  mine 
hunting.  The  high  resolution  necessary  for  classification 
of  mines,  however,  requires  a  very  long  array.  One  of  the 
most  attractive  ways  to  satisfy  these  competing  require¬ 
ments  is  through  aperture  synthesis. 

The  principle  of  aperture  synthesis  is  the  coherent 
combination  of  successive  returns  from  a  transmitter  lo¬ 
cated  on  a  moving  platform.  This  coherent  processing  of 
received  echoes  provides  enhanced  resolution  in  azimuth 
compared  with  standard  sonar  systems.  The  key  to  SAS  sig¬ 
nal  processing  is  estimating  the  time  delay  of  each  re¬ 
ceived  pulse  while  compensating  for  motion  of  the  array  and 
the  target.  This  process  can  be  thought  of  as  focusing  of 
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the  phase  history  of  the  signal  [1] .  Traditional  sonar 
arrays  use  multiple  transmit  and  receive  elements  to  simul¬ 
taneously  transmit  and  then  simultaneously  receive  acoustic 
signals.  Conversely,  synthetic  aperture  arrays  synthesize 
the  array  by  processing  the  output  electrical  signals  of  a 
single  transmit  and  receive  element  or  group  of  elements  as 
the  element  travels  through  the  medium.  A  comparison  of 
the  traditional  array  and  a  synthetic  aperture  can  be  seen 
in  Figure  1 . 


Traditional  Array  Synthetic  Aperture  Array 


Figure  1 .  Traditional  and  synthetic  aperture  arrays 
[After  Ref .  2 ] . 


SAS  systems  are  similar  to  Synthetic  Aperture  Radar 
(SAR)  systems  in  many  respects  but  SAS  signal  processing 
provides  some  unique  challenges.  SAR  systems  have  wider 
absolute  bandwidth  and  higher  carrier-to-bandwidth  ratio, 
or  quality  factor  Q  .  The  high  Q  nature  of  SAR  dictates 
that  typical  motion  and  range  estimation  errors  have  little 
effect  on  the  echo.  The  errors  result  in  a  shift  of  a  few 
centimeters  on  a  signal  envelope  of  a  few  meters.  In  con¬ 
trast,  SAS  motion  and  range  estimation  errors  of  a  few 
centimeters  result  in  a  shift  of  the  target  envelope  that 
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is  also  on  the  order  of  a  few  centimeters  [3] .  The  strong 
dependence  of  SAS  systems  on  accurate  range  and  motion 
estimation  is  the  focus  of  the  research  presented  in  this 
thesis . 

Several  other  challenges  must  be  considered  when  de¬ 
signing  a  SAS  system.  The  speed  of  the  platform  must  be 
high  enough  to  maintain  control  of  the  submerged  vehicle 
but  slow  enough  to  prevent  undersampling  of  the  target 
area.  As  described  above,  SAS  system  resolution  can  be 
significantly  degraded  by  residual  motion  errors  of  the 
transmit  and  receive  platforms. 

B.  THE  OCEAN  AS  A  LINEAR,  TIME -VARIANT,  SPACE-VARIANT 

FILTER 

Computer  simulation  of  acoustic  ray  propagation  in  an 
unbounded  ocean  medium  was  the  primary  means  of  evaluating 
the  exact  propagation  equations  and  comparing  these  tech¬ 
niques  to  common  approximations.  The  mathematical  frame¬ 
work  used  for  these  simulations  models  the  ocean  as  a  lin¬ 
ear,  time-variant,  space-variant  filter.  The  theory,  meth¬ 
ods,  and  derivation  of  this  linear  systems  approach  to 
acoustic  field  theory  are  contained  in  [4] .  A  block  dia¬ 
gram  model  of  pulse  propagation  is  shown  in  Figure  2.  The 
simulations  conducted  in  this  thesis  focus  on  determination 
of  exact  and  approximate  acoustic  propagation  paths  in¬ 
cluded  in  the  complex  frequency  response  of  the  ocean  me¬ 
dium.  The  complex  frequency  response  of  the  ocean  medium 
at  time  t ,  position  r,  and  frequency  /  due  to  application 

of  a  unit-amplitude  impulse  at  position  r0  is  HM  [t,r \fr0)  ,  as 
seen  in  Figure  2 . 
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Figure  2.  Block  diagram  model  of  small-amplitude  pulse 
propagation  [From  Ref.  5] . 

Models  used  for  simulation  are  based  on  bistatic  scat¬ 
tering  with  a  stationary  target.  This  scenario  is  intended 
to  model  the  performance  of  a  submerged  vehicle  using  a 
side  scan  sonar  to  image  stationary  mines.  The  transmitter 
and  receiver  are  co-located  on  the  platform.  Due  to  motion 
of  the  vehicle  the  position  of  the  transmitter  and  receiver 
is  time  varying.  It  is  important  to  note  that  the  vehicle 
will  be  in  motion  during  the  entire  transmit  and  receive 
cycle.  Some  SAS  signal-processing  algorithms  assume  that 
the  vehicle  is  stationary  during  each  transmit  and  receive 
cycle  and  moves  only  between  each  cycle  [6] .  This  simpli¬ 
fication  requires  motion  compensation  to  be  applied  during 
signal  processing  and  can  introduce  errors  in  the  image. 

The  models  used  in  this  thesis  are  exact  and  require  no 
motion  compensation.  The  exact  models  do,  however,  assume 
that  the  velocity  of  the  platform  is  constant  during  trans¬ 
mission  and  reception  of  individual  pulses. 

C .  RESEARCH  GOALS 

This  thesis  focuses  on  the  direct  comparison  of  exact 
and  approximate  signal  propagation  models  as  applied  to  a 
SAS  system  employed  to  image  stationary  targets.  This 
research  validates,  through  computer  simulation,  the  exact 
propagation  model.  The  exact  model  is  then  used  to  evalu- 
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ate  errors  introduced  by  approximate  methods.  The  goal  of 
this  research  is  to  provide  a  proven  SAS  signal-processing 
algorithm  that  removes  errors  due  to  propagation  range, 
angle,  and  platform  motion  estimation.  With  reduced  error, 
the  synthetic  aperture  array  has  increased  array  gain  and 
improved  image  quality. 


D .  DOD  RELEVANCE 

Mine  detection  and  classification  is  one  of  the  most 
important  challenges  faced  by  the  Navy.  Simple,  inexpen¬ 
sive,  effective  mines  are  available  to  all  our  potential 
adversaries.  Deployment  of  water-borne  mines  by  an  adver¬ 
sary  precludes  the  introduction  of  ships  or  troops  into  an 
otherwise  unprotected  area. 

The  effectiveness  of  mine  warfare  was  proven  in  the 
1991  Persian  Gulf  War.  The  USS  Tripoli  sustained  3.5  mil¬ 
lion  dollars  in  damage  after  being  struck  by  a  1500  dollar 
mine.  Later  in  the  war,  the  USS  Princeton  was  disabled  by 
a  similar  mine.  Repairs  to  the  USS  Princeton  cost  24  mil¬ 
lion  dollars  [7] .  Battle  group  commanders  are  justified  in 
their  reluctance  to  employ  assets  where  the  threat  of  wa¬ 
ter-borne  mines  exists. 

To  address  the  asymmetric  threat  of  mines,  the  U.S. 
Navy  is  currently  developing  a  Long-Term  Mine  Reconnais¬ 
sance  System  (LMRS) .  The  LMRS  is  an  unmanned  underwater 
vehicle  that  will  use  a  forward-looking  sonar  system  for 
mine  detection  and  a  side-scan  sonar  system  for  localiza¬ 
tion  and  classification.  The  side-scan  sonar  is  used  as 
the  transmitter  and  receiver  for  a  SAS  system.  This  system 
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provides  no  risk  access  to  mined  waters  while  collecting 
information  essential  to  mine  hunting  operations.  A  con¬ 
ceptual  drawing  of  the  LMRS  can  be  seen  in  Figure  3. 


Figure  3.  LMRS  conceptual  drawing  [From  Ref.  8] . 

E .  SUMMARY 

Reliable  detection  and  classif ication  of  water-borne 
mines  is  an  essential  capability  for  the  U.S.  Navy.  SAS 
systems  on  UUV' s  are  currently  under  development  to  address 
this  challenge.  Traditional  SAS  beamforming  algorithms  use 
approximations  to  calculate  acoustic  signal  propagation 
ranges  and  angles  of  incidence.  Exact  methods  of  determin¬ 
ing  these  parameters  have  been  developed  and  improve  system 
performance.  This  thesis  uses  computer  simulations  of 
exact  and  approximate  propagation  models  based  on  linear 
systems  theory  to  compare  the  performance  of  these  methods. 
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F.  THESIS  ORGANIZATION 

This  thesis  is  organized  into  five  chapters.  Chapter 
I  provides  an  overview  of  the  SAS  concept,  military  rele¬ 
vance  of  this  research,  and  a  general  description  of  the 
modeling  methods  and  evaluation  techniques  used  for  system 
simulation.  An  overview  of  the  application  of  linear  sys¬ 
tems  theory  to  acoustic  propagation  modeling  is  presented 
in  Chapter  II.  Chapter  III  introduces  the  mathematical 
models  used  for  simulation  of  SAS  as  applied  to  the  imaging 
of  stationary  water-borne  targets.  The  results  of  computer 
simulations  designed  to  compare  different  methods  of  esti¬ 
mating  acoustic  signal  propagation  ranges  and  received 
phase  angles  are  presented  in  Chapter  IV.  Finally,  Chapter 
V  presents  the  thesis  conclusions  and  recommendations  for 
future  research  in  SAS  signal  processing. 


The  next  chapter  describes  the  treatment  of  the  ocean 
as  a  linear,  time-variant,  space-variant  filter.  The  SAS 
mine-hunting  geometry  and  scattering  model  are  also  intro¬ 
duced  . 
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II.  THE  OCEAN  AS  A  LINEAR,  TIME -VARIANT ,  SPACE- 

VARIANT  FILTER 


This  chapter  presents  an  overview  of  the  application 
of  linear  systems  theory  to  acoustic  propagation  modeling. 


A.  INTRODUCTION 

The  computer  simulations  presented  in  this  thesis  are 
based  on  a  linear  systems  model  of  the  ocean  environment. 

A  linear  systems  approach  to  acoustic  pulse  propagation  is 
presented  in  [4]  and  [5] .  All  simulation  models  in  this 
thesis  are  based  on  treatment  of  the  ocean  as  a  linear 
filter . 

Propagation  of  small  amplitude  acoustic  pulses  in  the 
ocean  can  be  described  by  a  linear  wave  equation.  This 
linear  wave  equation  accurately  describes  the  propagation 
of  acoustic  pulses  from  the  source  to  the  discrete  point 
scatterer  and  from  the  discrete  point  scatterer  to  the 
receiver.  To  model  this  propagation,  we  treat  the  ocean  as 
a  linear,  time-variant,  space-variant  filter. 


B.  COMPLEX  FREQUENCY  RESPONSE  OF  THE  OCEAN 

The  propagation  of  small  amplitude  acoustic  pulses 
through  the  ocean  medium  can  be  described  by  the  linear, 
three-dimensional,  inhomogeneous  wave  equation,  [5] 

1  d 2 

v2>v(f>r) — 2,  V  3,2  yM(t’r)  =  xii(t,r)  ,  (2.1) 

c  (r)  at 
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where  xM(t,r)  is  the  source  distribution  at  time  t  and  posi¬ 
tion  r  with  units  of  inverse  seconds  and  vM(hr)is  the  veloc¬ 
ity  potential  of  the  acoustic  field  at  time  t  and  position 
r  with  units  of  squared  meters  per  second.  The  vector 
r=(x,y,z)  describes  the  location  of  any  point  in  three- 
dimensional  space.  The  position  dependent  speed  of  sound 
in  meters  per  second  is  described  by  c(r)  .  Note  that  both 
the  source  distribution  and  the  velocity  potential  are 
functions  of  time  and  position.  The  time-variant  property 
of  the  linear  filter  allows  us  to  model  motion  of  the 
transmitter,  target,  and  receiver.  This  ability  to  accu¬ 
rately  model  motion  of  all  platforms  is  essential  in  SAS 
modeling  where  transmitter  and  receiver  motion  is  used  to 
synthesize  the  array.  This  time-variant  property  also  is 
used  to  account  for  changes  in  the  ocean  medium  with  time. 
The  space-variant  property  allows  for  the  presence  of 
boundaries,  discrete  point  scatterers  and  spatial  variation 
in  ambient  density  and  the  speed  of  sound. 

The  solution  of  Eq.  (2.1)  can  be  obtained  by  treating 
the  ocean  as  a  linear,  time-variant,  space-variant  filter. 
Figure  4  is  the  linear  system  block  diagram  of  the  ocean  as 
a  time-variant,  space-variant  filter.  The  input-output 
relationship  is  given  by  [5] 

y M  r)  =  L  ^  ’  r«  )kM  (*>  '■ r  \f0’  r0  )dt0  dr0  r  (2.2) 

where  hM  (/,r  | f 0 , r0 )  is  the  time-variant,  space-variant,  impulse 
response  (Green's  function)  of  the  system.  The  impulse 
response  hM[t,r\t0,r^  is  the  response  of  the  filter  at  time  t 
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and  position  r  =  (x,y,z)  due  to  application  of  a  unit  ampli¬ 
tude  impulse  at  time  t0  and  position  r0  =  (x0,_y0,z0)  .  The  input 
acoustic  signal  xM(t, r)  is  the  source  distribution  at  time  t 
and  position  r  and  yM(t, r)  is  output  acoustic  signal  at  time 
t  and  position  r.  It  is  important  to  realize  that  Eq. (2.2) 
is  a  four-fold  integral  since  dr0=  dx0  dy0  dz0  . 


xM(t,  r) 


>'w(Ar) 


Figure  4.  A  linear,  time-variant,  space-variant,  ocean 
medium  filter  [From  Ref.  5] . 


Figure  4  also  shows  the  time-variant,  space-variant, 
transfer  function  of  the  ocean  filter  HM  (?,r  |/,v)  defined  by 
[5] 


|/,r0)exp[+y2^v*(r-r0)]Jr0,  (2.3) 

where  /  represents  input  frequency  components  in  Hertz  and 
V  =  (/Y,/7,/z)  is  a  vector  representation  of  the  spatial  fre¬ 
quencies  in  cycles  per  meter  in  the  X,  Y,  and  Z  directions, 
respectively.  The  time-variant,  space-variant,  complex 
frequency  response  of  the  ocean  at  frequency  /  is  [5] 

Hu{t,r  |/,r0)  =  J  \t0,r0)exip[-j27rf(t-t0)\dt0  .  (2.4) 
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The  strategy  followed  in  Reference  5  is  to  first  find 
the  time  harmonic  solution  to  the  linear,  inhomogeneous, 
three-dimensional  wave  equation.  The  time  harmonic  solu¬ 
tion  is  then  used  to  find  a  pulse  solution  to  the  wave 
equation  by  using  Fourier  transform  techniques. 

The  focus  of  this  thesis  is  on  the  complex  frequency 
response  of  the  ocean  medium,  HM  (hr  |/,r0)  ,  as  defined  by 

Eq.  (2.4)  .  The  complex  frequency  response  of  the  ocean 
medium  integrates  the  characteristics  of  the  propagation 
path  between  transmitter,  discrete  point  scatterer,  and 
receiver.  Simulations  were  conducted  using  both  exact  and 
approximate  propagation  ranges  to  account  for  signal  propa¬ 
gation  delay  as  determined  by  propagation  path  of  the 
acoustic  field. 

C .  SUMMARY 

The  ocean  medium  can  be  modeled  as  a  linear,  time- 
variant,  space-variant,  filter.  This  approach  allows  us  to 
develop  a  time-variant,  space-variant,  complex  frequency 
response  of  the  ocean  that  accurately  describes  acoustic 
signal  propagation  and  therefore  provides  a  solution  to  the 
linear,  inhomogeneous,  wave  equation  given  by  equation 
(2.1) .  In  the  next  chapter,  the  complex  frequency  response 
of  the  ocean  is  derived  for  the  specific  case  of  SAS  ap¬ 
plied  to  stationary  object  detection. 
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III.  MODELING  OF  BISTATIC  SCATTERING  FOR  SYNTHETIC 

APERTURE  SONAR 


This  chapter  introduces  the  mathematical  model  used 
for  simulation  of  synthetic  aperture  sonar  (SAS)  as  applied 
to  the  imaging  of  stationary  water-borne  targets.  Several 
traditional  approximations  and  an  exact  solution  for  propa¬ 
gation  ranges  are  introduced  and  the  limitations  of  each 
are  discussed. 

A.  INTRODUCTION 

The  main  advantage  of  a  SAS  system  is  the  ability  to 
synthesize  an  extremely  long  array  while  maintaining  a 
relatively  small  actual  size  of  the  physical  array.  Syn¬ 
thetic  aperture  sonar  beamforming  and  signal  processing 
relies  on  properly  steering  and  focusing  the  aperture  beam 
pattern  in  order  to  co-phase  all  the  received  signals.  Due 
to  the  effects  of  motion  in  the  SAS  problem,  the  propaga¬ 
tion  path  from  the  transmitter  to  the  discrete  point  scat- 
terer  and  from  discrete  point  scatterer  to  the  receiver  is 
time  varying.  Traditionally,  simple  approximations  are 
used  to  determine  these  propagation  ranges  and  angles  of 
incidence  and  scatter.  This  chapter  introduces  the  exact 
bistatic  scattering  model  for  SAS  systems  presented  in  [9] 
as  well  as  some  common  approximations  used  when  modeling 
acoustic  signal  propagation  in  SAS  systems. 
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B. 


BISTATIC  SCATTERING  MODEL 


1.  Assumptions 

The  focus  of  this  section  is  to  develop  the  complex 
frequency  response  of  the  ocean  for  a  SAS  system  with  sta¬ 
tionary  targets.  In  the  binomial  approximation  models  and 
the  exact  propagation  models  described  in  this  thesis,  the 
speed  of  sound  and  the  ambient  density  of  the  ocean  are 
treated  as  constants.  This  allows  us  to  treat  the  ocean 
medium  as  homogeneous.  As  a  result,  sound  rays  will  travel 
in  straight  lines  from  the  transmitter  to  the  discrete 
point  scatterer  and  from  the  discrete  point  scatterer  to 
the  receiver.  Propagation  models  involving  surface  and 
bottom  interactions  can  be  developed  in  an  identical  manner 
to  those  presented  for  propagation  between  the  transmitter 
and  discrete  point  scatterer  and  discrete  point  scatterer 
and  receiver  [9] . 

The  transmitter  and  receiver  velocity  vectors  are  as¬ 
sumed  to  be  constant  only  during  the  transmission  and  re¬ 
ception  of  the  acoustic  pulse.  The  platform  may  change  its 
speed  at  any  time  in  any  direction  between  transmission  and 
reception  and  between  cycles.  This  feature  eliminates  the 
need  for  motion  compensation  typically  applied  to  most  SAS 
signal  processing  algorithms.  The  binomial  approximation 
and  exact  propagation  models  assume  only  that  the  position 
and  velocity  vectors  of  the  transmitter  and  receiver  can  be 
determined  during  signal  transmission  and  reception.  De¬ 
termination  of  the  vehicle  velocity  vector  is  not  a  trivial 
problem,  however,  and  accurate  determination  of  vehicle 
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position  is  an  area  of  extensive  research  in  SAS  signal 
processing  and  vehicle  design  [1,  10,  11]  . 

In  order  to  allow  comparison  of  the  various  methods  of 
propagation  range  determination,  the  scattering  function  of 

the  discrete  point  scatterer,  gl  (/, n0l, n"2 ) ,  is  assumed  to  be 

unity . 


2 .  Propagation  Geometry 

We  begin  our  development  of  the  complex  frequency  re¬ 
sponse  of  the  ocean  by  defining  the  velocity  vectors  for 
the  transmitter,  V0 ,  the  discrete  point  scatterer,  V, ,  and 
the  receiver,  V2  .  These  velocity  vectors  are  given  by  [9] 

\  =  V0nVor  (3.1) 

\=V,  hVi,  (3.2) 

and 

V2=F2>V2,  (3.3) 

where  V0 ,  Vl ,  and  V2  are  the  speeds  in  meters  per  second  of 
the  transmitter,  discrete  point  scatterer,  and  receiver, 
respectively.  The  dimensionless  unit  vectors  nr  ,  nr  ,  and 

77, v  define  the  directions  of  V0 ,  V,,  and  V2 ,  respectively. 

Velocity  vectors  given  by  Eqs .  (3.1)  through  (3.3)  are  con¬ 
stant  as  discussed  in  Section  III.B.l.  Motion  is  consid¬ 
ered  to  begin  at  time  t  =  tm. 

The  time-varying  geometry  of  the  scattering  model  is 
modified  using  relative  velocity  vectors  to  allow  the 
transmitter  and  the  discrete  point  scatterer  to  be  treated 
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as  being  motionless.  The  velocity  vector  of  the  discrete 
point  scatterer  relative  to  the  velocity  vector  of  the 
transmitter  in  the  direction  of  the  velocity  vector  of  the 
discrete  point  scatterer  hv  ,  is  [9] 

Vi.o  —  Vj  —  (%[  *  ~  V\~  Vo  [nVx  •  hv^  jj  / hVi  .  (3.4) 

The  initial  propagation  model  allows  for  motion  of  the 
discrete  point  scatterer,  so  we  define  a  relative  velocity 
vector,  V21  .  The  vector  V21  is  the  velocity  vector  of  the 

receiver  relative  to  the  velocity  vector  of  the  discrete 
point  scatterer  in  the  direction  of  the  velocity  vector  of 
the  receiver  hr  and  is  given  by  [9] 

V2,i  —  V2  —  (nv^  •  V,  j  hv^  —  V2  —  Vl  (hy2  •  hy  nVi  .  (3.5) 

The  position  vectors  from  the  origin  to  the  transmit¬ 
ter,  discrete  point  scatterer,  and  receiver  when  motion 
begins  at  t  =  tm  are  given  by  r0  =  (x0,y0,z0) ,  r \=(xl,yl,zl) ,  and 

r2  =(x2,y2,z2) ,  respectively.  Using  these  position  vectors,  we 

can  further  define  the  position  vector  from  the  transmitter 
to  the  discrete  point  scatterer  as  [9] 

ro,i  =  ri  - ro  (3-6) 

and  the  position  vector  from  the  discrete  point  scatterer 
to  the  receiver  as  [9] 

ri,2  =  r2  —  ri  -  (3.7) 

The  dimensionless  unit  vectors  from  source  to  discrete 
point  scatterer  and  from  discrete  point  scatterer  to  re¬ 
ceiver  are  [9] 
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and 


‘0,1 


0,1 


‘0,1 


(3.8) 


‘1,2 


(3.9) 


respectively.  The  ranges  in  meters  from  the  transmitter  to 
discrete  point  scatterer  and  from  the  discrete  point  scat- 
terer  to  receiver  are  given  by  [9] 

ro,\  ~  |ro,i  |  (3.10) 

and 

rU2  =|ru| ,  (3.11) 


respectively.  This  bistatic  scattering  geometry  is  shown 
in  Figure  5 . 

It  is  now  important  to  define  the  parameters  of  our 
model  with  respect  to  the  time  instance  when  the  acoustic 
field  is  first  incident  upon  the  discrete  point  scatterer 
and  the  time  instance  when  the  acoustic  field  is  first 
incident  upon  the  receiver.  The  transmitted  acoustic  field 
is  first  incident  upon  the  discrete  point  scatterer  at  some 
time  t'  seconds  where  t'>tm  .  This  allows  us  to  define  the 
time  difference  from  transmission  of  the  acoustic  field  to 
first  incidence  upon  the  discrete  point  scatterer  as  [9] 

A  t'=t'-tm,  (3.12) 


Similarly,  the  time  instant  when  the  scattered  acoustic 
field  is  first  incident  upon  the  receiver  is  t>t'>tm.  The 
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time  difference  between  t'  and  first  incidence  of  the  acous¬ 
tic  field  upon  the  receiver  is  given  by  [9] 


A t"  =  t-t\  t>  t’>t  .  (3.13) 

~  m  '  ' 


Figure  5.  Bistatic  scattering  geometry  when  motion  begins 
at  time  t  =  tm  seconds.  Point  0,  i^(r0),  is  the  transmitter; 

point  1,  is  the  discrete  point  scatterer;  and  point 
2,  P2(r 2),  is  the  receiver.  All  three  platforms  are  in  mo¬ 
tion  [From  Ref.  5] . 


By  referring  to  Figure  6  we  can  express  the  position 
vector  from  the  point  source  to  the  discrete  point  scat¬ 
terer  at  time  t'  as  and  the  position  vector  from  the 
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discrete  point  scatterer  to  the  receiver  at  time  t  as  r"2  . 

The  position  vector  from  the  origin  to  the  discrete  point 

/  ,  /  , 

scatterer  at  time  t  is  ^  .  The  position  vector  from  the 

rr 

origin  to  the  receiver  at  time  t  is  r2  . 

3 .  Velocity  Potential  and  Complex  Frequency  Response 

Our  goal  is  to  develop  an  expression  for  the  complex 
frequency  response  of  the  ocean  in  order  to  compare  differ¬ 
ent  methods  for  estimation  of  propagation  range  from  trans¬ 
mitter  to  discrete  point  scatterer  and  from  discrete  point 
scatterer  to  receiver.  Prediction  of  the  propagation 

ranges  |r^ 1 1  and  |r*2|  will  allow  us  to  correctly  steer  and 

focus  the  beam  pattern  of  our  synthetic  array  resulting  in 
the  co-phasing  of  the  output  electrical  signals  at  the 
different  receiver  locations.  Once  a  model  for  the  complex 
frequency  response  is  developed,  we  can  predict  the  acous¬ 
tic  field  (velocity  potential)  at  the  receiver. 

When  the  transmitted  acoustic  field  is  first  incident 
upon  the  discrete  point  scatterer  at  time  t' ,  the  position 
vector  from  the  transmitter  to  the  discrete  point  scatterer 
is  given  by  [9] 

ro,i=ro,i+A^  Vo  •  (3.14) 


19 


So<5(r  -  r,' ) 


Figure  6.  Bistatic  scattering  geometry  when  the  transmit¬ 
ted  acoustic  field  is  first  incident  upon  the  discrete 
point  scatterer  at  time  t'  seconds  and  when  the  scattered 
acoustic  field  is  first  incident  upon  the  receiver  at  time  t 
seconds  where  t>t'>tm.  Point  0,  i^,(r0),  is  the  transmitter; 

point  1,  is  the  discrete  point  scatterer,  and  point 

2,  P2[ r,),  is  the  receiver.  All  three  platforms  are  in  mo¬ 

tion  [From  Ref.  5] . 
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Recall  that  we  are  modeling  the  propagation  of  the  acoustic 
field  through  the  ocean  as  transmission  through  a  linear 
time-variant,  space-invariant  filter.  Therefore,  the 
acoustic  field  incident  upon  the  discrete  point  scatterer 

/  /  /  /  /  /\ 
at  time  t  and  position  =  I xx  ,y1  ,zx  I  is  [9] 

Jiu  (*'>  )  =  S()Hm  (t\  r'u  \f )  exp(+j2xft') ,  t'>tm  (3.15) 


where 


tx 


0,1 


exp| 

-jk  i 

r0,i  ) 

4  K 

r0,i 

(3.16) 


is  the  time-variant,  space-invariant,  complex  frequency 
response  of  the  ocean  at  frequency  /  hertz.  The  source 
strength  S0  is  given  in  cubic  meters  per  second  and  k  is 
the  wavenumber  in  radians  per  meter  given  by 

k  =  27Tf/c  =  In IX  .  (3.17) 


An  identical  method  is  used  in  [9]  to  develop  an  ex¬ 
pression  for  the  acoustic  field  incident  upon  the  receiver 
at  time  t  .  The  position  vector  between  the  discrete  point 
scatterer  and  the  receiver  at  time  t>t'>tm  is  given  by  [9] 


<2=ru+A/(V2 -V1>0)  +  AfXi 


(3.18) 


The  acoustic  field  (velocity  potential)  incident  upon  the 
receiver  at  time  t>t'>t  is  [9] 

m  L  J 

yM  ri  )  =  s»hm  (k  r*  \f,  r0 )  exp  {+j2nft ) ,  t>t>  tm  (3.19) 


where 
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The  complex  scattering  function  of  the  discrete  point 
scatterer,  gl  {f,n01,fi"2} ,  is  dependent  upon  frequency  and  upon 

the  angles  of  incidence  and  scatter.  For  simplicity,  this 
function  is  assumed  to  be  unity  in  this  thesis  so  that  we 
may  independently  investigate  the  phase  term  in  Eg. (3.20) . 


C.  PROPAGATION  RANGE  DETERMINATION 

The  key  to  proper  co-phasing  of  the  output  electrical 
signals  from  the  SAS  array  is  estimation  of  the  received 
phase  of  the  acoustic  field  incident  upon  the  receiver. 
Examination  of  Eq. (3.20)  reveals  that  the  phase  term, 

exp|^-/'£(|r(,)1|  +  |r"2|)  ,  depends  directly  on  frequency  and  the 

magnitude  of  the  position  vectors  and  r*2  .  It  is  in  the 

solution  for  these  time-varying  ranges  that  this  thesis  is 
focused.  Accurate  determination  of  the  propagation  range 
of  the  acoustic  signal  is  especially  important  in  recent 
SAS  projects  such  as  the  Long-Term  Mine  Reconnaissance 
System  (LMRS)  where  relatively  high  carrier  frequencies 
will  be  used  to  increase  image  resolution.  Small  errors  in 
propagation  range  estimation  can  introduce  large  errors  in 
signal  phase  estimation  at  the  receiver,  significantly 
degrading  the  array  performance. 
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1 .  The  Stop-and-Hop  Approximation  and  the  Moving  Re¬ 
ceiver  Correction 

One  approach  to  estimation  of  propagation  ranges  is  to 
assume  that  the  transmitter  and  receiver  are  stationary 
during  signal  transmission  and  reception.  This  assumption 
is  acceptable  for  many  synthetic  aperture  radar  (SAR)  sys¬ 
tems,  but  due  to  the  slow  propagation  speed  of  acoustic 
waves  in  water  relative  to  electromagnetic  waves  in  air, 
this  assumption  is  invalid  for  a  SAS  system  imaging  targets 
at  any  significant  range  [6] .  An  initial  correction  may  be 
applied  by  compensating  for  the  motion  of  the  receiver  in 
the  time  interval  between  signal  transmission  and  signal 
reception.  A  method  for  received  phase  correction  due  to  a 
moving  receiver  is  presented  in  [6] .  Figure  7  shows  the 
geometry  involved  in  this  correction.  The  initial  range 
from  transmitter  to  target  is  R{  .  The  initial  range  from 

target  to  receiver  is  R2  .  The  distance  traveled  by  the 
receiver  during  signal  propagation  is  Vt'  where  V  is  the 
vehicle  speed  and  t'  is  the  total  propagation  time  between 
transmission  and  reception  of  the  acoustic  signal.  The 
propagation  time  t'  is  an  estimation  of  the  total  propaga¬ 
tion  time  from  transmitter  to  target  and  from  target  to 
receiver  based  on  a  reference  range  to  the  scene  center. 

The  range  i?3  describes  the  range  from  target  to  receiver  at 

time  t'  .  Each  received  acoustic  signal  is  multiplied  by  the 
phase  correction  term  [6] 


exp 


+ j—AR 

c 


(3.21) 
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Figure  7.  The  moving  receiver  correction  for  the  stop- 
and-hop  geometry  [After  Ref.  6] . 

where 

A  R  =  R3-R2.  (3.22) 

Compensating  for  the  motion  of  the  receiver  during 
signal  propagation  does  remove  some  error  inherent  in  the 
stop-and-hop  assumption  but  several  additional  factors  are 
not  considered.  The  motion  of  the  transmitter  and  receiver 
during  signal  transmission  and  reception  is  ignored.  In 
addition,  the  vehicle  velocity  vector  is  assumed  to  be 
constant  for  the  entire  duration  of  signal  propagation. 

This  correction  further  assumes  that  the  total  propagation 
time,  t'  in  Figure  7,  from  transmitter  to  target  and  target 
to  receiver  can  be  approximated  by  calculating  the  propaga¬ 
tion  time  to  some  reference  range,  R0  ,  such  that  [6] 

,  2  R. 

t  = — .  (3.23) 

c 
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The  reference  range,  R0  ,  is  the  range  to  the  center  of 
the  image  area  of  the  SAS  system  [6] .  Throughout  the  re¬ 
mainder  of  this  thesis,  the  moving  receiver  correction  is 
applied  to  the  stop-and-hop  model  and  referred  to  as  the 
stop-and-hop  approximation  model. 

We  can  define  the  stop-and-hop  model  presented  in  [6] 
in  terms  of  the  bistatic  scattering  model  presented  in 
Section  III.B  of  this  thesis.  The  range  from  transmitter 
to  discrete  point  scatterer  is  given  by 

Rx  - 11',  —  r0|  ,  (3.24) 

where  rt  and  r0  are  the  position  vectors  from  the  coordinate 
system  origin  to  the  discrete  point  scatterer  and  transmit¬ 
ter,  respectively  (see  Fig.  5) .  The  time-varying  range 
from  discrete  point  scatterer  to  receiver  is  given  by 

(f)  =  |r2(f)  —  r,|,  t>  t  ,  (3.25) 

where  r2  (^)  is  the  time-varying  range  from  the  coordinate 
system  origin  to  the  receiver  given  by 

r2(t)  =  r2(t')  +  \t,  t>  t' .  (3.26) 

The  vector  V  is  the  velocity  vector  of  the  receiver,  and  t' 
is  given  by  Eq.  (3.23)  . 

2 .  The  Binomial  Approximation  Method 

The  magnitude  of  the  position  vectors  and  x"X2  can  be 

approximated  by  using  the  binomial  expansion  presented  in 
[5] .  The  magnitude  of  the  position  vector  from  the  trans- 
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mitter  to  the  discrete  point  scatterer  at  time  t'  can  be 
rewritten  as  [5] 


<1  =r0,VT+& 


f  b  b  ^ 

1  + - +••• 


v 


2  8 


b  <1 


J 


(3.27) 


where 


b  =  2 


Vytf 


0,1 


(«o,i  *  hv) 


+  - 


2  r, 


0,1 


(3.28) 


Recall  that  in  our  bistatic  scattering  model,  t'  is  the  time 
instant  when  the  transmitted  acoustic  field  is  first  inci¬ 
dent  upon  the  discrete  point  scatterer  and  that  t'>tm.  For 
the  case  of  a  SAS  system  imaging  a  stationary  target. 


V,  =0 


(3.29) 


Therefore,  V10  given  by  Eg. (3.4)  is  undefined.  As  a  result. 


1,0 


=  0 


(3.30) 


and 


vlt0  =r1>0  =  o 


(3.31) 


In  this  case,  b  given  by  Eq.  (3.28)  is  equal  to  zero  and 
Eq. (3.27)  reduces  to 


(3.32) 


which  is  exact. 


The  binomial  approximation  can  also  be  applied  to  the 
estimation  of  |r*2|  .  The  range  |r*2|  can  be  written  as  [5] 
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1,2  1,2 


'Jl  +  b 


( 


'1,2 


,  b  b2 

1  + - +  ••• 

2  8 


b\<\ 


(3 . 33) 


where 


b  =  —  {[»U  •  (■ V,  -  v, )]  A •  V, , )  A,'}  . 

ri,2 


(3.34) 


The  expression  for  b  can  be  simplified  using 
Eqs .  (3.29)  and  (3.30)  for  the  stationary  target.  Thus, 
Eg.  (3.5)  reduces  to 


V  =  V 
’  2,1  ’2 


(3 . 35) 


and  Eg.  (3.34)  can  be  rewritten  as 


6  =  — {(A/  +  A^)(»12*V2)}  . 


(3.36) 


1,2 


Using  Eqs.  (3.12)  and  (3.13), 

At'  +  At"  =  t'-tm  + 1  - 1'  =  t  -  tm  =  A t 

The  expression  for  b  now  reduces  to 


(3.37) 


b  =  —  {A^(«i,2#V2)} 


(3.38) 


1,2 


Equation  (3.33)  is  only  valid  if  |Z>|<1.  Note  that 

k2.V2|<|V2|.  (3.39) 

Now,  if  we  impose  the  more  stringent  criterion  that 

2|VJ  Af 

^ — <0.1,  (3.40) 

then  the  binomial  approximation  can  be  simplified  by  using 
only  the  first  two  terms  of  Eq.  (3.33)  .  Doing  so  yields 
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(3.41) 


‘1,2 


!  \l  + 


‘1,2 


'(V2-V10)  A/'+(n12.V2l)A/'. 


Equation  (3.41)  can  be  further  simplified  using  Eqs . (3.30), 
(3.35),  and  (3.37)  such  that 

|rn2 1  =  '1,2  +  («1,2  •  v2 )  Af  .  (3.42) 


3 .  Exact  Propagation  Range  Method 

Reference  9  presents  a  method  for  determining  the 

time-varying  propagation  ranges  |rj j (?)|  and  |r*2(V)|  exactly  for 

the  general  case  of  bistatic  scattering  when  all  three 
platforms  (transmitter,  target,  and  receiver)  are  in  motion 
and  for  the  special  case  of  a  SAS  system  imaging  a  station¬ 
ary  target.  The  constant  value  of  range  from  transmitter 
to  discrete  point  scatterer  when  the  transmitted  acoustic 
field  is  first  incident  upon  the  discrete  point  scatterer 
at  time  instant  t'  >tm  is  given  by  [9] 


1 0,1 


'0,1 


(3.43) 


The  constant  value  of  range  between  discrete  point 
terer  and  the  receiver  when  the  scattered  field  is 
incident  upon  the  receiver  at  time  instant 
by  [9] 


scat- 
first 
is  given 


K  ±  A3  +4AQ 

2A0 


(3.44) 


where 


4  = 


(3.45) 
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(3.46) 


and 


B,=2r 


1,2 


fh>+2,.  d 

0,1  2  1 

c  c 


n  —  „2  ^0  , 

Co  —  ro,l  2  +  ^  r0,l  ri,2 


+  r, 


1,2 


(3.47) 


The  time-varying  propagation  ranges  are  given  by  [9] 

|  ru(0 1  =  [ri,2  +  2  ru  («i,2  *  % ) A t  +  V02  (At)2  ] 1/_ ,  t>tm  +  T,  (3.48) 

and 

|ro.i(0|  =  cA/-|r'2(0|,  t>tm+T,  (3.49) 

where 


c  c 


is  the  time  delay  in  seconds  or  the  amount  of  time  it  takes 
for  the  transmitted  acoustic  signal  to  begin  to  appear  at 
the  receiver  after  motion  begins  at  time  instant  tm,  and 


A  t  =  t-t 


t  >  t.„  +  T 


(3.51) 


D .  SUMMARY 

Accurate  estimation  of  the  propagation  range  from  the 
transmitter  to  target  and  from  target  to  receiver  is  essen¬ 
tial  to  the  proper  co-phasing  of  the  output  electrical 
signals  from  the  different  receiver  locations.  The  tradi¬ 
tional  stop-and-hop  method,  corrected  for  receiver  motion, 
can  be  used  to  estimate  propagation  ranges.  Alternately, 
the  propagation  ranges  can  be  determined  using  the  bistatic 
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scattering  model  presented  in  this  chapter.  The  propaga¬ 
tion  ranges  |rj j(?)|  and  |r^2(0|  can  be  estimated  using  a  bino¬ 
mial  approximation  or  calculated  exactly.  The  next  chapter 
presents  the  results  of  computer  simulations  used  to  com¬ 
pare  the  accuracy  of  the  stop-and-hop  approximation,  bino¬ 
mial  approximation,  and  the  exact  propagation  range  calcu¬ 
lation  methods. 
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IV.  PROPAGATION  RANGE  AND  RECEIVED  PHASE  ANGLE 
ESTIMATION  COMPARISON 


This  chapter  presents  the  results  of  computer  simula¬ 
tions  designed  to  compare  different  methods  of  estimating 
acoustic  signal  propagation  ranges  and  received  phase  an¬ 
gles  for  a  SAS  system  imaging  a  stationary  target.  Simula¬ 
tion  results  of  the  receiver  output  electrical  signals  for 
a  SAS  array  with  Fast  Fourier  Transform  (FFT)  beamforming 
and  various  propagation  range  estimation  methods  are  also 
presented . 


A.  INTRODUCTION 

1 .  Method  of  Evaluation 

Accurate  estimation  of  the  propagation  range  from  the 
transmitter  to  the  discrete  point  scatterer  and  from  the 
discrete  point  scatterer  to  the  receiver  is  essential  to 
properly  co-phase  the  output  electrical  signals  in  the  SAS 
system.  Using  the  models  presented  in  Chapter  III,  the 
simulations  presented  in  this  chapter  compare  the  propaga¬ 
tion  ranges  and  received  phase  angles  computed  using  the 
binomial  approximation  method  and  stop-and-hop  approxima¬ 
tion  method  to  the  exact  propagation  ranges  and  exact  re¬ 
ceived  phase  angles.  Comparisons  are  based  on  the  percent 
error  in  range  from  transmitter  to  discrete  point  scat¬ 
terer,  percent  error  in  range  from  discrete  point  scatterer 
to  receiver,  and  percent  error  in  the  estimation  of  the 
phase  of  the  signal  incident  upon  the  receiver.  All  simu¬ 
lations  were  conducted  using  MATLAB  version  6.0. 
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Comparisons  of  propagation  range  and  received  phase 
estimation  methods  are  performed  for  a  single  trans¬ 
mit/receive  cycle.  The  propagation  range  error  for  the 
time-varying  range  from  transmitter  to  discrete  point  scat- 
terer  is  given  by 


Exact 


<i(0 


Approx 


(4.1) 


Similarly,  the  propagation  range  error  for  the  time-varying 
range  from  discrete  point  scatterer  to  receiver  is  given  by 

|r'(0|  -|r'2(0|  •  (4.2 

I  1’2V  '\ Exact  I  J'2V  '\ Approx 


The  percent  range  error  for  the  time-varying  range  from 
transmitter  to  discrete  point  scatterer  is  given  by 


<i(0 

Exact 

<i(0 

Approx 

<i(0 

Exact 

xl00% 


(4.3) 


The  percent  range  error  for  the  time-varying  range  between 
discrete  point  scatterer  and  receiver  is  given  by 

kr2(o  L 

— - ^Lxl00%.  (4.4) 

Exact 


Finally,  the  phase  in  radians  of  the  acoustic  signal  inci¬ 
dent  upon  the  receiver  is  given  by 


and  the  percent  phase  error  is  given  by 

e(t)  -e(t) 

v  /  Exact  v  '  Approx 


9(t) 


Xl00% 


Exact 


(4.5) 


(4.6) 
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2. 


Simulation  Parameters 


The  computer  simulations  developed  for  this  thesis 
simulate  vehicle  motion  in  any  direction  in  an  unbounded, 
homogeneous,  ocean  medium.  Although  the  simulation  places 
no  restrictions  on  target  location,  the  target  is  always 
assumed  to  be  at  some  depth  below  the  transmit /receive 
platform.  Unless  otherwise  noted,  the  vehicle  begins  mo¬ 
tion  at  time  tm  at  the  origin  of  the  coordinate  system  and 

is  imaging  a  target  forward  of  the  beam  (0  =  45°)  at  a  range, 
R  (see  Fig.  8) .  For  simplicity,  platform  motion  is  along 
the  positive  Z  axis. 


Z 


Figure  8.  SAS  simulation  geometry  [After  Ref.  5] . 
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Propagation  ranges  are  calculated  with  the  equations 
presented  in  III.C.l  through  III.C.3  using  the  assumptions 
presented  in  III.B.l. 

The  platform  speeds  and  target  ranges  used  for  com¬ 
parison  are  similar  to  those  that  have  been  specified  for 
the  Navy's  LMRS  system  [2] .  The  LMRS  side-scan  sonar  sys¬ 
tem  is  used  to  synthesize  the  SAS  array.  The  side-scan 
sonar  array  covers  both  sides  of  the  six-foot  long  LMRS 
platform,  allowing  independent  imaging  on  both  the  port  and 
starboard  sides  of  the  vehicle.  The  vehicle  is  designed  to 
travel  between  3  and  7  knots  [2] .  The  minimum  speed  is 
based  on  maintaining  maneuverability.  The  maximum  speed  is 
set  to  ensure  adequate  sampling  of  the  target  area.  Figure 
9  shows  the  relationship  between  physical  array  length, 
platform  speed,  and  maximum  range. 

MAXIMUM  SAS  COVERAGE  (PER  SIDE)FOR  VARIOUS 


PLATFORMSPEEDS  AMD  ARRAY  LENGTHS 


Figure  9.  LMRS  Range  and  Speed  limitations  [From  Ref.  2] . 
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3. 


FFT  Beamforming 


The  simulations  presented  in  the  next  section  compare 
the  accuracy  of  the  approximate  stop-and-hop  propagation 
range  method  and  the  approximate  binomial  approximation 
propagation  range  method  for  each  transmit /receive  cycle  of 
the  platform.  To  provide  further  comparison  and  to  better 
analyze  the  effects  of  approximation  errors  on  SAS  systems, 
the  propagation  range  calculations  from  the  stop-and-hop 
approximation,  binomial  approximation,  and  exact  methods 
are  used  to  estimate  the  phase  weights  in  a  Fast  Fourier 
Transform  (FFT)  beamforming  algorithm. 

To  independently  investigate  the  performance  of  each 
method  due  to  phase  errors,  we  assume  the  scattering  func¬ 
tion  of  the  target  is  unity  as  discussed  in  Sections 
III.B.l  and  III.B.3.  We  also  neglect  the  propagation  range 

term,  \6n2  |r(^a  |  |r^2 1 ,  in  the  denominator  of  Eq.(3.20).  Neglect¬ 
ing  this  amplitude  factor  allows  us  to  focus  on  the  effect 
of  propagation  range  error  on  received  phase  prediction. 

Three  simple  sinusoids,  each  at  a  different  frequency, 
are  used  to  simulate  acoustic  propagation.  As  a  result, 

Eq.  (3.19)  can  be  expressed  as 


yM  (*>r2'(0)  =  XexP [-1%  (|<i(0|  +  |<2(0|)]  Qxp{j2nft ) , 

1=1 


(4.7) 


where 


c 


is  the  wavenumber  for  each  different  frequency  and  fl  =  10 
kHz,  f2  =  20  kHz,  and  /3  =  30  kHz. 


35 


We  assume  that  the  output  electrical  signal  of  the  re¬ 
ceiver  is  directly  proportional  to  the  acoustic  field  inci¬ 
dent  upon  the  receiver.  Therefore, 

=  .  (4.9) 

We  now  describe  the  FFT  beamforming  algorithm  used  in 
this  thesis  and  is  shown  in  Fig.  10.  FFT  beamforming  is 
accomplished  by  first  taking  the  time-domain  Fourier  trans¬ 
form  of  the  output  electrical  signals  at  the  receiver  at 
different  positions  as  the  SAS  system  images  a  target  at  a 
known  location.  These  signals  are  simulated  using  the 
exact  propagation  model  to  compute  the  acoustic  signal 

incident  upon  the  array  given  by  Eq.  (4.7),  where  |r^ j (f )|  and 

Ir^  (?)|  are  given  by  Eqs .  (3.49)  and  (3.48),  respectively.  The 

resulting  output  electrical  signals  are  given  by  Eq. (4.9) . 
These  exact  electrical  signals  are  then  phase  weighted  with 
phase  weights  computed  using  one  of  the  approximate  meth¬ 
ods  . 

These  approximate  phase  weights  are  computed  by  first 
determining  the  acoustic  field  incident  upon  the  receiver 

using  Eq.  (4.7),  where  the  values  for  |%(?)|  and  |r1/'2(f)|  are 

computed  using  either  the  binomial  approximation  method  or 
the  stop-and  hop  approximation  method.  The  approximate 

stop-and-hop  values  for  |r(' : (f )|  and  |r,//2(/)|  are  given  by 
Eqs.  (3.24)  and  (3.25),  respectively.  The  approximate  bino¬ 
mial  approximation  values  for  |rj j(?)|  and  |rf2(y)|  are  given  by 

Eqs.  (3.32)  and  (3.42),  respectively.  The  negative  of  the 
phase  spectrum  of  the  time-domain  Fourier  transform  of  the 
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T/R  Cycle  N  T/R  Cycle  2  T/R  Cycle 


output  electrical  signals  given  by  Eq.  (4.9)  determine  the 
appropriate  approximate  phase  weights  to  use  to  phase 
weight  the  receiver  output  electrical  signals. 


Figure  10.  FFT  Beamforming  Algorithm. 

The  approximate  phase  weights  are  then  multiplied  by 
the  time-domain  Fourier  transform  of  the  output  electrical 
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signals  determined  using  the  exact  method.  The  exact  out¬ 
put  electrical  signals  have  now  been  phase  weighted  using 
approximate  phase  weights.  We  then  compute  the  inverse 
Fourier  transform  of  the  phase-weighted  signals.  The  re¬ 
sulting  time-domain  signals  are  added  to  produce  the  total 
output  electrical  signal  from  the  SAS  array  given  by 

yT  (t)  =  X  ^C„  (?  ~(tmn  +T„))'  (4.10) 

n= 1 

where  N  is  the  total  number  of  transmit /receive  cycles  used 
to  create  the  SAS  array  and  yc  [t  —  {tm  +  T(i)j  is  the  output 

electrical  signal  from  the  receiver  for  the  n‘h  transmit/ 
receive  cycle  after  phase  weighting.  The  signals  are  syn¬ 
chronized  in  time  by  subtracting  the  time  delay  for  each 

transmit /receive  cycle,  [tm  +  Tn),  where  tm  is  the  time  at  the 

beginning  of  the  rih  cycle  and  Tn  is  given  by  Eq.  (3.50)  .  In 
the  FFT  beamforming  algorithm  shown  in  Fig.  10,  we  assume 
that  the  output  electrical  signals  have  been  synchronized 
in  time  before  performing  the  forward  FFT.  With  the  trans¬ 
mit/receive  platform  constantly  in  motion,  the  initial 
range  from  the  transmitter  to  the  discrete  point  scatterer 
and  from  the  scatterer  to  the  receiver  will  be  different  at 
the  beginning  of  each  transmit /receive  cycle.  These  ini¬ 
tial  ranges  are  updated  to  reflect  the  geometry  at  the 
beginning  of  each  cycle  and,  as  a  result,  the  propagation 
delay,  Tn  ,  is  recomputed  for  each  transmit /receive  cycle. 

The  FFT  beamforming  procedure  is  described  in  detail  in 
[4]  . 
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B. 


PROPAGATION  RANGE  AND  RECEIVED  PHASE  COMPARISON 


1 .  Stop-and-Hop  With  Moving  Receiver  Correction 

The  performance  of  the  stop-and-hop  propagation  model 

is  evaluated  using  Eqs .  (4.1)  through  (4.6)  where  r'(7)  is 

I  I  Exact 

given  by  Eq.  (3.4  9),  t (^)|^  is  given  by  Eq.  (3.24),  t 

is  given  by  Eq.  (3.48)  ,  and  is  given  by  Eq.  (3.25)  . 

Figure  11  shows  the  errors  introduced  by  the  stop-and- 
hop  approximation  for  a  vehicle  speed,  V0 ,  of  2-7  knots, 
i?  =  213  yards,  and  frequency,  /,  of  100  kHz  for  a  single 
transmit /receive  cycle.  In  this  case,  the  stop-and-hop 
approximation  with  the  moving  receiver  correction  intro¬ 
duces  no  error  in  the  estimation  of  the  range  from  discrete 

point  scatterer  to  receiver,  |ip,(h)|  •  The  error  in  the  esti¬ 
mation  of  range  from  transmitter  to  discrete  point  scat¬ 
terer,  |rjj(/)|  ,  is,  however,  significant.  The  error  intro¬ 
duced  by  the  stop-and-hop  method  results  in  a  maximum  per¬ 
cent  received  phase  error  of  -0.5%  at  F0=7  knots. 

At  the  minimum  platform  speed  of  3  knots,  we  can  in¬ 
vestigate  estimation  errors  from  the  stop-and-hop  approxi¬ 
mation  over  all  possible  target  ranges  observing  the  limi¬ 
tations  described  in  Figure  9.  Figure  12  shows  propagation 
range  and  received  phase  error  for  F0=3  knots, 

/  =  100  kHz,  with  R  varied  from  25-450  yards.  From  Figure 
12  we  can  see  that  the  estimation  error  for  |r^  1(?)|  and  the 
subsequent  phase  error  increases  as  the  range  to  the  target 
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decreases.  This  suggests  that  the  stop-and-hop  approxima¬ 
tion  provides  unacceptable  performance  inside  some  minimum 
range . 


|r’0  1  (t)|  Error:  Stop-and-Hop 


lr’o  iWI  %  Error:  Stop-and-Hop 


x1g4  Phase  Error:  Stop-and-Hop 


%  Phase  Error:  Stop-and-Hop 


Figure  11.  Propagation  range  and  received  phase  error 
using  the  stop-and-hop  approximation  for  F0=2-7  knots, 
/  =  100  kHz,  and  i?  =  213  yards. 
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|r'Q  1(t)|  Error:  Stop-and-Hop  |r'0 1  (t)|  %  Error:  Stop-and-Hop 


X104  Phase  Error:  Stop-and-Hop  %  Phase  Error:  Stop-and-Hop 


Figure  12.  Propagation  range  and  received  phase  error 
using  the  stop-and-hop  approximation  for  V0=3  knots, 

/  =  100  kHz,  and  i?  =  25-450  yards. 
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It  is  interesting  to  examine  the  accuracy  of  the  stop- 
and-hop  approximation  with  the  target  at  different  initial 
angles  with  respect  to  the  array.  The  angle  (j)  is  a  cylin¬ 
drical  coordinate  measured  from  array  broadside  to  the 
target  projection  in  the  X-Z  plane  at  the  time  motion  be¬ 
gins,  as  seen  in  Figure  8.  This  simulation  is  conducted 
such  that  R  remains  constant  at  300  yards  while  (/)  is  var¬ 
ied  from  0°  to  90°  .  The  target  depth  is  set  at  150  yards 
below  the  platform.  The  platform  speed,  V0 ,  is  5  knots. 
Figure  13  shows  that  propagation  range  and  received  phase 
error  increase  as  (f)  increases.  This  is  an  important  obser¬ 
vation  since  all  SAS  systems  employ  the  coherent  addition 
of  multi-aspect  data  to  form  high-resolution  images.  The 
ability  to  accurately  image  a  target  at  higher  aspect  an¬ 
gles  results  in  improved  along-track  image  resolution  [2] . 
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x  1  q4  Phase  Error:  Stop-and-Hop  %  Phase  Error:  Stop-and-Hop 


Figure  13.  Propagation  range  and  received  phase  error 
using  the  stop-and-hop  approximation.  Angle  <j)  varies  from 

0°  to  90°.  Speed  V0  =  5  knots,  /  =  100  kHz,  and  i?  =  300  yards. 
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2 .  Binomial  Approximation  Method 

The  performance  of  the  binomial  approximation  method 
for  propagation  range  estimation  is  evaluated  using  Eqs . 

(4.1)  through  (4.6)  where  r'(7)  is  given  by  Eq.(3.49), 
1^(7)  is  given  by  Eq.  (3.32),  ^(Y)  is  given  by  Eq. 

I  '  I  Approx  I  1  I  Exact 

(3.48),  and  |r'/,(7)|  is  given  by  Eq.(3.42). 

Simulations  examining  the  accuracy  of  the  binomial  ap¬ 
proximation  model  were  conducted  in  the  same  manner  as 
those  for  the  approximate  stop-and-hop  model.  Once  again, 
the  binomial  approximation  estimates  of  the  propagation 

ranges  |rj j(?)|  and  and  the  resulting  estimated  phase  of 

the  acoustic  signal  incident  upon  the  receiver  were  com¬ 
pared  with  the  exact  solutions  for  these  parameters. 

Figure  14  shows  the  errors  introduced  by  the  binomial 
approximation  for  a  vehicle  speed,  V0 ,  of  2-7  knots, 

R  =  213  yards,  and  frequency,  /,  of  100  kHz  for  a  single 
transmit /receive  cycle.  In  this  case,  the  binomial  ap¬ 
proximation  introduces  extremely  small  error  in  the  estima¬ 
tion  of  the  range  from  the  transmitter  to  the  discrete 

point  scatterer,  |rj j (?)|  .  Maximum  estimation  errors  for  the 
binomial  approximation  of  |rj j(?)|  are  on  the  order  of  10-6 
meters  at  V0=7  knots.  The  error  from  the  estimation  of 
range  from  discrete  point  scatterer  to  receiver,  |r*2(7)|,  is 
also  extremely  small.  Maximum  estimation  errors  for  the 
binomial  approximation  of  |r*2(0|  are  on  the  order  of  1CT2 
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meters  at  V0  =7  knots.  The  resulting  phase  error  introduced 
by  the  binomial  approximation  method  results  in  a  maximum 
percent  received  phase  error  of  0.004%  at  F0=7  knots. 
Recalling  that  the  stop-and-hop  approximation  resulted  in  a 
maximum  percent  phase  error  of  -0.5%  for  the  same  scenario, 
it  is  obvious  that  the  binomial  approximation  greatly  re¬ 
duces  propagation  range  and  received  phase  estimation  er¬ 
ror  . 

At  the  minimum  platform  speed  of  3  knots,  we  can  again 
investigate  estimation  errors  from  the  binomial  approxima¬ 
tion  over  all  possible  target  ranges  observing  the  limita¬ 
tions  described  in  Figure  9.  Figure  15  shows  propagation 
range  and  received  phase  error  for  V0=3  knots  and  /  =  100 
kHz  while  R  is  varied  from  25-450  yards.  From  Figure  15  we 
can  see  that  the  estimation  error  for  |r^ j (f)|  increases  as  the 
target  range,  R  ,  increases.  Conversely,  the  range  estima¬ 
tion  error  for  |r*2(/)|  decreases  as  R  increases,  and  is  much 

greater  than  the  range  estimation  error  for  |rj x (?)|  in  all 

cases.  This  results  in  a  phase  error  that  increases  as  R 
decreases,  a  result  similar  to  that  observed  for  the  stop- 
and-hop  approximation.  This  suggests  that  the  binomial 
approximation  also  may  provide  unacceptable  performance 
inside  some  minimum  range.  Comparing  the  stop-and-hop 
approximation  method  with  the  binomial  approximation 
method,  we  again  see  that  the  binomial  approximation  method 
introduces  much  less  error  in  the  estimation  of  propagation 
ranges  and  received  phase. 
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Figure  14.  Propagation  range  and  received  phase  error 
using  the  binomial  approximation  method  for  V0  =  2  —  7  knots, 
/  =  100  kHz,  and  i?  =  213  yards. 
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Figure  15.  Propagation  range  and  received  phase  error 
using  the  binomial  approximation  method  for  V0=3  knots, 
/  =  100  kHz,  and  i?  =  25-450  yards. 
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We  can  also  examine  the  accuracy  of  the  binomial  ap¬ 
proximation  method  with  the  target  at  different  initial 
angles  with  respect  to  the  array.  The  angle  (j)  is  a  cylin¬ 
drical  coordinate  measured  from  array  broadside  to  the 
target  projection  in  the  X-Z  plane,  as  seen  in  Figure  8. 
This  simulation  is  conducted  such  that  R  remains  constant 
at  300  yards  while  </)  is  varied  from  0"  to  90°  .  The  target 
depth  is  set  at  150  yards  below  the  platform.  The  platform 
speed,  V0=5  knots.  Figure  16  shows  that  propagation  range 
and  received  phase  error  decrease  as  (j)  increases.  This 
suggests  that  the  binomial  approximation  method  would  be 
able  to  more  accurately  combine  the  multi-aspect  data  in¬ 
herent  to  SAS  systems.  Once  again,  we  can  see  that  the 
binomial  approximation  introduces  a  maximum  of  0.002%  phase 
error  compared  to  -0.25%  phase  estimation  error  per  trans¬ 
mit/receive  cycle  introduced  by  the  stop-and-hop  approxima¬ 
tion. 
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Figure  16.  Propagation  range  and 
using  the  binomial  approximation, 
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C.  FFT  BEAMFORMING  SIMULATIONS 

To  provide  some  insight  into  the  performance  of  the 
various  propagation  models  when  applied  to  a  SAS  system, 
each  model  was  used  to  calculate  the  phase  weights  for  use 
in  an  FFT  beamforming  algorithm  as  described  in  IV. A. 3, 
where  fx  — 10  kHz,  /2  =  20  kHz,  and  /3=30  kHz. 

The  array  is  synthesized  as  the  vehicle  travels  at  the 
designated  speed  in  a  straight  line  along  the  positive  Z 
axis  . 

Figure  17  shows  results  of  FFT  beamforming  with 
(a)  F0=3  knots  and  (b)  F0=7  knots.  The  array  is  synthe¬ 

sized  using  10  transmit /receive  cycles  with  the  target  at 
an  initial  range,  R  ,  of  213  yards  and  ^  =  45°.  Phase  weight¬ 
ing  the  output  electrical  signals  using  the  binomial  ap¬ 
proximation  method  to  determine  received  phase  of  the  inci¬ 
dent  acoustic  signal  at  both  3  and  7  knots  produces  a  wave¬ 
form  that  is  almost  exact.  The  errors  in  the  stop-and-hop 
approximation,  however,  are  significant  and  phase  weighting 
the  output  electrical  signals  using  this  method  results  in 
very  low  array  gain.  The  poor  performance  of  the  stop-and- 
hop  approximation  results  from  the  relatively  high  error  in 
the  received  acoustic  signal  phase  estimate. 
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Synthetic  Aperture  Sonar  Total  Output  Electrical  Signal 


Synthetic  Aperture  Sonar  Total  Output  Electrical  Signal 


Figure  17.  FFT  beamforming  with 
V0  =  7  knots.  The  array  is  synthes 
receive  cycles  with  the  target  at 
213  yards  and  0  =  45°. 


(a)  V0=3  knots  and  (b) 

ized  using  10  transmit/ 
an  initial  range,  R  ,  of 
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Synthetic  Aperture  Sonar  Total  Output  Electrical  Signal 


Figure  18.  FFT  beamforming  with  F0=3  knots  and 
(a)  i?  =  100  yards,  (b)  R  =  250  yards,  and  (c)  R  =  500  yards. 
The  initial  target  angle,  0  =  45°.  The  array  is  synthesized 
using  10  transmit /receive  cycles. 
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Figure  18  shows  the  results  of  FFT  beamforming  with 
V0=3  knots  and  (a)  i?  =  100  yards,  (b)  i?  =  250  yards,  and  (c) 

i?  =  500  yards.  The  initial  target  angle,  0  =  45°  in  all  three 
cases  and  the  array  is  synthesized  using  10  trans¬ 
mit/receive  cycles.  We  again  see  that  phase  weighting  the 
output  electrical  signals  using  the  binomial  approximation 
of  received  acoustic  signal  phase  is  almost  exact.  Phase 
weighting  using  the  stop-and-hop  approximation  method  for 
estimating  received  acoustic  signal  phase  results  in  very 
low  array  gain  in  all  cases. 

The  errors  introduced  by  the  binomial  approximation 
method  seen  in  Figures  14,  15,  and  16  are  not  large  enough 
to  significantly  affect  the  steering  and  focusing  of  the 
beam  pattern  for  typical  vehicle  speeds  and  target  ranges. 
When  the  array  is  placed  within  1  yard  of  the  target  and 
the  vehicle  is  allowed  to  travel  at  7  knots,  however,  the 
binomial  approximation  errors  become  evident.  Figure  19 
shows  the  results  of  FFT  beamforming  with  R= 1  yard, 

V0=7  knots,  and  0  =  45°  for  an  array  synthesized  with  10 
transmit /receive  cycles.  In  this  case,  the  errors  intro¬ 
duced  in  the  estimation  of  propagation  ranges  with  the 
binomial  approximation  method  are  a  significant  percentage 
of  the  total  range. 
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Synthetic  Aperture  Sonar  Total  Output  Electrical  Signal 


Figure  19.  FFT  beamforming  with  R  =  1  yard,  VQ  =  1  knots,  and 

0  =  45°  for  an  array  synthesized  with  10  transmit/ 
receive  cycles. 
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D .  SUMMARY 

The  propagation  ranges  and  received  acoustic  signal 
phase  can  be  estimated  using  both  the  stop-and-hop  method 
and  the  binomial  approximation  method.  Both  methods  intro¬ 
duce  some  error  into  the  calculation  of  propagation  ranges 
and  received  acoustic  signal  phase.  Simulation  results 
indicate  that  the  error  introduced  by  the  binomial  approxi¬ 
mation  is  significantly  lower  than  the  error  introduced  by 
the  stop-and-hop  approximation  for  typical  LMRS  vehicle 
speeds  and  target  ranges.  Estimates  of  received  phase  from 
both  the  stop-and-hop  approximation  and  the  binomial  ap¬ 
proximation  methods  were  used  to  phase  weight  the  output 
electrical  signals  in  an  FFT  beamforming  algorithm.  The 
errors  introduced  by  the  stop-and-hop  approximation  method 
significantly  degrade  the  performance  of  the  FFT  beamform¬ 
ing  algorithm.  Conversely,  the  binomial  approximation 
method  for  received  acoustic  signal  phase  estimation  can  be 
used  to  accurately  phase  weight  the  output  electrical  sig¬ 
nals  at  typical  LMRS  vehicle  speeds  and  target  ranges.  The 
thesis  conclusions  are  presented  in  the  next  chapter. 
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V. 


SUMMARY  AND  CONCLUSIONS 


This  chapter  provides  a  summary  of  the  propagation 
range  and  received  acoustic  signal  phase  estimation  meth¬ 
ods.  Based  on  the  simulations  presented  in  the  last  chap¬ 
ter,  the  accuracy  of  the  stop-and-hop  approximation  and  the 
binomial  approximation  methods  is  analyzed.  Finally,  fur¬ 
ther  research  in  SAS  acoustic  signal  processing  is  recom¬ 
mended  . 

A.  SUMMARY 

Accurate  estimation  of  acoustic  signal  propagation 
range  and  received  phase  is  critical  to  SAS  system  perform¬ 
ance.  To  determine  the  phase  of  the  acoustic  signal  inci¬ 
dent  upon  the  receiver  at  all  locations  along  the  synthetic 
array,  the  time-varying  range  from  the  transmitter  to  the 
target  and  from  the  target  to  the  receiver  must  be  deter¬ 
mined.  The  estimated  phase  of  the  received  acoustic  signal 
can  then  be  used  to  steer  and  focus  the  beam  pattern  of  the 
synthetic  aperture  array. 

The  SAS  concept  is  currently  being  applied  to  the 
problem  of  water-borne  mine  classification  and  imaging  in 
the  Navy's  Long-Term  Mine  Reconnaissance  System  (LMRS) . 

This  requires  high  resolution  imaging  at  long  ranges. 

Under  these  conditions,  small  errors  in  the  estimates  for 
propagation  range  and  subsequent  phase  estimates  can  sig¬ 
nificantly  degrade  the  SAS  system's  ability  to  steer  and 
focus  the  array  beam  pattern. 

This  thesis  analyzed  the  performance  of  the  stop-and- 

hop  approximation  and  the  binomial  approximation  models  for 
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determining  propagation  range  and  received  phase  of  acous¬ 
tic  signals  for  a  SAS  system  imaging  a  stationary  target. 
The  stop-and-hop  method,  presented  in  [6],  initially  as¬ 
sumes  that  the  vehicle  is  stationary  during  transmission, 
propagation,  and  reception  of  the  acoustic  signal.  A  sim¬ 
ple  correction  is  applied  to  account  for  receiver  motion. 
The  binomial  approximation  method,  presented  in  [5],  devel¬ 
ops  equations  for  the  range  from  transmitter  to  target  and 
from  target  to  receiver  for  the  general  problem  of  bistatic 
scattering.  Approximate  solutions  for  propagation  range 
equations  are  obtained  by  performing  a  binomial  expansion 
of  the  propagation  range  equations.  The  derivations  of 
this  technique  place  some  restrictions  on  the  simulation 
parameters,  but  at  the  same  time  ensure  minimal  estimation 
error . 

The  stop-and-hop  approximation  method  and  the  binomial 
approximation  method  were  compared  using  simulations  that 
exactly  determine  the  acoustic  signal  propagation  range  and 
received  phase  for  a  SAS  system  imaging  a  stationary  tar¬ 
get.  This  exact  method,  presented  in  [9],  for  the  first 
time  allows  us  to  investigate  the  accuracy  of  common  ap¬ 
proximations  and  their  effect  on  SAS  system  performance. 

This  thesis  also  presented  an  overview  of  the  model 
used  to  simulate  the  propagation  of  acoustic  signals  for  a 
SAS  system.  In  the  simulations  presented,  the  ocean  is 
treated  as  a  linear,  time-variant,  space-variant  filter. 

The  complex  frequency  response  of  the  ocean  filter  is  the 
primary  focus  of  the  simulations  presented. 
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B.  CONCLUSIONS 

Simulations  presented  in  this  thesis  are  based  on 
typical  LMRS  vehicle  speeds  and  target  ranges.  Performance 
of  the  different  approximation  methods  was  compared  for  a 
single  transmit /receive  cycle  and  then  the  estimates  of 
received  acoustic  signal  phase  were  used  to  compare  the 
performance  of  each  method  applied  to  an  FFT  beamforming 
algorithm . 

The  stop-and-hop  approximation  introduced  a  maximum  of 
-0.5%  received  phase  error  with  the  vehicle  traveling  at 
the  LMRS  maximum  specified  speed  and  at  the  maximum  target 
range.  The  error  in  propagation  range  and  received  phase 
estimation  increased  as  target  range  increased,  resulting 
in  a  maximum  of  -2.5%  received  phase  estimation  error  with 
the  transmit /receive  platform  at  the  minimum  LMRS  speed  and 
the  target  at  25  yards.  The  estimation  errors  introduced 
by  the  stop-and-hop  approximation  method  are  small,  but 
significant.  This  is  best  demonstrated  when  the  stop-and- 
hop  approximation  method  is  used  to  determine  the  phase 
weights  in  an  FFT  beamforming  algorithm  applied  to  a  SAS 
array.  The  error  introduced  by  the  stop-and-hop  approxima¬ 
tion  method  significantly  degraded  the  performance  of  the 
FFT  beamforming  algorithm.  It  was  also  evident  that  the 
accuracy  of  the  stop-and-hop  approximation  method  decreased 
as  the  target  aspect  angle  increased.  This  suggests  that 
the  stop-and-hop  approximation  method  is  poorly  suited  to 
modern  SAS  techniques  that  combine  multi-aspect  data  to 
improve  image  resolution. 

The  binomial  approximation  method  estimated  acoustic 
signal  propagation  range  and  received  phase  per  trans- 
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mit/receive  cycle  with  much  less  error  than  the  stop-and- 
hop  approximation  method.  At  the  maximum  LMRS  speed  and 
target  range,  the  received  phase  estimation  error  was  ap¬ 
proximately  100  times  less  than  the  error  introduced  by  the 
stop-and-hop  approximation.  Similarly,  at  low  platform 
speeds  and  close  target  ranges  the  estimation  error  intro¬ 
duced  by  the  binomial  approximation  method  was  1000  times 
less  than  the  error  introduced  by  the  stop-and-hop  method. 
In  fact,  the  errors  introduced  by  the  binomial  approxima¬ 
tion  were  insignificant  when  the  estimates  of  received 
acoustic  signal  phase  were  used  to  phase  weight  the  output 
electrical  signals  of  the  array  in  an  FFT  beamforming  algo¬ 
rithm.  Only  at  very  close  target  ranges  (less  than  1  yard) 
and  high  platform  speeds  did  the  errors  introduced  by  the 
binomial  approximation  significantly  affect  the  performance 
of  the  SAS  system.  In  contrast  to  the  stop-and-hop  ap¬ 
proximation  method,  the  accuracy  of  the  binomial  approxima¬ 
tion  method  improved  as  the  target  aspect  angle  increased. 

The  propagation  range  and  received  acoustic  signal 
phase  can  be  determined  exactly  for  the  SAS  system  using 
the  techniques  presented  in  [9] .  These  exact  solutions  do 
not  add  significant  complexity  to  the  SAS  signal-processing 
problem  while  eliminating  errors  introduced  by  approximate 
methods.  These  exact  solutions  place  no  restrictions  on 
target  range  and  vehicle  velocity,  allowing  for  maximum 
flexibility  as  SAS  systems  evolve. 
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C.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

This  thesis  considered  improvements  to  SAS  system 
beamforming  using  the  exact  and  binomial  approximation 
methods  for  determining  propagation  range  and  received 
acoustic  signal  phase  versus  using  the  stop-and-hop  ap¬ 
proximation.  The  methods  presented  in  [9]  also  provide 
exact  solutions  for  the  time-varying  angles  of  incidence 
and  scatter  at  the  discrete  point  scatterer  and  the  time- 
varying  angles  of  incidence  at  the  receiver.  This  addi¬ 
tional  information  should  be  incorporated  into  models  of 
the  scattering  functions  for  mines  and  mine-like  objects 
since  scattering  functions  are  not  only  functions  of  fre¬ 
quency,  but  also  functions  of  the  angles  of  incidence  and 
scatter.  Scattering  functions  should  be  included  in  future 
simulations  in  order  to  determine  their  impact  on  SAS  sys¬ 
tem  performance.  Using  LMRS  operational  parameters,  typi¬ 
cal  oceanographic  data,  and  real  world  target  simulations, 
and  end-to-end  SAS  simulation  should  be  developed.  The 
simulations  would  be  critical  to  the  development  of  SAS 
operational  guidelines,  capabilities,  and  future  improve¬ 
ment  . 
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